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Abstract 

Intermittency of £^(x,g) = exp[<7|S , 7v(x)| 2 ] as N — > +oo is investigated on a d-dimensional 
torus A, when Sn(x) is a finite Steinhaus series of (2iV + l) d terms normalized to (\Sn(x)\ 2 ) = 1. 
Assuming ergodicity of £^{x,g) as N — ► +oo in the domain g < 1, where limN^ +00 (£]y(g)) 
exists, transition to intermittency is proved as g increases past the threshold gth = 1- This tran- 
sition goes together with a transition from (assumed) ergodicity at g < gth to a regime where 
liin/v^+ooIIAK^A^g))] -1 J a £n{x, g) d d x = at g > gth- In this asymptotic sense one can say that 
ergodicity is lost as g increases past the value g = 1. 
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I. INTRODUCTION 



This paper is the first of a series devoted to studying intermittency of the solution to the 
random PDE, 

d t £ N (x,t) - 2^A£jv(:r,t) = X\S N (x,t)\ 2 £ N (x,t), 
t>0, x e A C R d , and £ N (x, 0) = 1, 

as N — > +00, where Sisr(x,t) is a sum of (2iV + l) d modes with i.i.d. random phases. Here 
A > is the coupling constant and m 7^ is a complex mass with Im(m) > 0. 

For Im(m) = and Re(m) 7^ 0, (1) models the scattering of an incoherent laser by an 
optically active medium. In this context, a simpler version in which Sn(x, t) is approximated 
by a Gaussian random field was first considered by Akhmanov et al. in nonlinear optics [1], 
and by Rose and DuBois in laser-plasma interaction [2]. The latter investigated the diver- 
gence of the average solution to (1) heuristically and numerically. The same problem was 
analyzed from a more rigorous mathematical point of view in [3], [4], and [5]. Going beyond 
the Gaussian approximation requires further specification of Sw(x,t). In realistic models 
laser light is represented by a superposition of a finite number of monochromatic beamlets 
with i.i.d. random phases [6]. The class of Sn(x, t) considered in (1) is a straighforward 
generalization of those models. For every N < +00, Sjy(x,t) is bounded and according to 
[5] there is no divergence of the average of £^{x ) t). To get interesting results from [5] one 
needs to work out the A" — > +00 limit of their theory. A possible alternative approach is 
suggested in the introduction of [3]. It is explained there that the divergence of the average 
solution to (1) indicates a change in the nature of £x{x,t) which undergoes a transition to 
intermittency. Taking it the other way round leads to characterize £n(x, t) by its intermit- 
tency, the divergence of its average taking a back seat. That is the approach followed in this 
work. Note that the problem (1) is complementary to the one considered in [5] in which A" 
is fixed and Sjy(x,t) is a sum of A" independent Gaussian random variables. By the central 
limit theorem, the results of [5] are expected to coincide with those of the present work in 
the limit A" — > +00. 

For a given t > 0, £n(x, t) is said to be intermittent if, for every integer p > 1, the space 
average of \£n{x, t)\ p is almost surely determined by higher and higher and more and more 
widely spaced peaks of 1^(^,01 as N — > +00. Intermittency of £^(x,t) can be inferred 
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1/p 



(2) 



from the almost sure chain of strong inequalities [7], 

K i^i / Mx,t)|d d a;< •••< t^t / IM^T^x 

I A I ./A LI A I A 

where /(iV) <C g(iV) means liminfAr^ +00 g(N)/f(N) = +oo. To prove that (2) does provide 
a sufficient condition for intermittency, choose for every realization for which (2) is fulfilled 
a sequence {f p (N)} such that 

' 1 ' * ' (3) 



i « Mm « im^oi d d x, 



« / P (iV) « 



|^(x,t)| p+1 ^x 



1 V(p+i) 



(4) 



and, for every integer p > 1, 

from (3) and (4) it follows that, Vp > 1, 

jT \E N {x, t)\ p l leN{x , t)l < fp _ l{N) d d x < U-^NY <^±.Jj£ N ( x , t )\Pd d x, 

hence 

W\ L ^ N<yX, ^ ' P dd% ~ W iv ^ l P1 i^( a; ' t )i >/ f- i ( 7v ) 

almost surely as TV — > +oo. Equations (5) means that when A" — > +oo, the space average 
of \£n(x, t)\ p is almost surely determined by the region of A in which |£jv(x, t)| > f p -i(N). 
Now, we must prove that this region gets smaller and smaller as A" — > +oo. This is easily 
done for p > 2 as it follows immediately from (4) that 



(5) 



(6) 



almost surely as A" — > +oo. For p = I assume that there exists e > small enough such 
that lim supjv^+oo |A| _1 f A \£n(x, t)\ e d d x < +oo with probability one, (this will be the case 
for the type of £n considered in this paper). Then, 
If, ,w , 1 



| A , l\ £N{x ,t)\ >M N)d d x < /o(iy 1 )£|A| \£n(x, t)\ £ d d x -> 0, 



(7) 



almost surely as A" — > +oo, which completes the proof that £jy(x,t) is intermittent when 
(2) holds. 

In this paper we determine the intermittency threshold defined as the smallest value of A 
above which £ N (x,t) is intermittent, assuming ergodicity as A" — > +oo for A small enough. 
By ergodicity as A" — > +oo we mean that if limjv_,. +00 (|£jv|) exists, then 



lim -4t / \£n(%, t)\ d d x = lim (\£n\), 

|A| A V 71 JV-H-oo 



(8) 



almost surely, where (•) denotes average over the realizations of Sn- We consider the simplest 
case |m| = +00, i.e. (1) without the ASn term, for a time-independent driver SN(x,t) = 
S N (x). 

The outline of the paper is as follows. In Section II we specify our model. Section III 
deals with the asymptotic behavior of (£n(x, t)) for large N. Intermittency of S^[x,t) and 
loss of ergodicity are investigated in Section IV. 

II. MODEL AND DEFINITIONS 

We assume that A is a d- dimensional torus of length L and volume |A| = L d . In the 
following we take L — 1 without loss of generality. For any given N e N, let An = 
{n e Z d : n G [-N, N] d }, Carder = (2iV + l) d , and assume that S N is of the form 

Sn ( x ) = (2N + iy/2 Yl exp + 27171 ' x ^ ' ( 9 ) 

^ ' neA N 

where the 9 n are i.i.d. random phases uniformly distributed over [0, 2tt[. In the terminology 
of the theory of random series of functions, Sn is called a Steinhaus (finite) series [8]. The 
average over the realizations of the 9 n is denoted by (•)$, or simply by (•) if there is no risk 
of confusion. 

In the limit \m\ = +00 and for Sn given by (9), the solution to (1) reduces to £isr(x,t) = 
exp(At|S , Ar(x)| 2 ). Introducing the average gain factor g = Xt(\S N (x)\ 2 ) = \t and using the 
fact that £n(%, t) is actually a function of x and g only, one is led to study the intermittency 
of the field 

£n(x,q) = exp [g\S N (x)\ 2 ] , (10) 

as N — > +00. The onset of intermittency will be characterized by the intermittency thresh- 
old, g t h, defined by 

g t h = inf{g > : £N(x,g) is intermittent}. (11) 

III. ASYMPTOTIC BEHAVIOR OF (£ N (x,g)) FOR LARGE N 

As we will see in the following, the intermittency properties of £N(x,g) depend on the 
behavior of (£n(x, <?)) for large N. This behavior is summarized in the following two lemmas. 



4 



Lemma 1 If g < 1, then V N > 0, (£n(x, g)) < (1 — g) 1 , and 

lim (S N (x,g)} = -^—. (12) 

1 — g 

Proof. Let h be a complex-valued zero-mean Gaussian random variable with (/i 2 ) = and 
{\h\ 2 } = (2N + l)~ d . Write exp [g\S N (x)\ 2 ] as 

exp [g\S N (x)\ 2 ] = (e^ 2N+l "> d ' 2 [SN{x)h* + s N { x yh\^ (13) 
Let u = \h\ 2 and define 

/» = «-lnIo(2 v ^), (14) 

where I is the modified Bessel function of zero order. From (9) and the integral represen- 
tation of I [9], one gets, 

(S N (x,g)) = ^ e V9^) d/2 lSN(x)^+s N ( x rh]^^ 

= ^ e V9(™+l) d / 2 [S N (x)h*+S N (x)*h]^ ^ 

= (27V + 1) " / + °° p e-W'^M 

^ J — oo J — oo 

r+oo 

= (2N + l) d / e -( 2JV+1 > d ^ u W (15) 
Jo 

If p < 1, it follows from (15) and the inequality In l Q (2y/gu) < gu that (£n(x, g)) is bounded 
above by 

(^(x, <?)> < (2N + l) d / e- {2N+1)d{1 - 9)u du = — — . (16) 

Jo i — g 

Furthermore, it can be easily checked that f g (u) is minimum at the boundary u — with 
/ 9 (0) = and f' g (0) = 1 — g > 0. The asymptotic behavior of (15) in the large N limit is 
thus determined by the vicinity of u = and one finds 

<m-,«)> = t ^[ 1 -°(^)_ 

hence (12), which completes the proof of Lemma 1. □ 



(AT^+oo), (17) 



Lemma 2 If g > I, then 3<y g > such that (£ N (x,g)) behaves like (2N + l) d / 2 exp[7 9 (2iV + 
l) d ] as N -> +oo. 

Proof It can be checked from (14) that if g > 1 there exists a unique number -u > such 
that / s (m) reaches its minimum at u — u , with f g (u ) < 0, f' g (u ) = 0, and f g (u ) > 0. 



Write 7 9 = —f g (u ) > 0. The asymptotic behavior of (15) in the large N limit is now 
determined by the vicinity of u — uo, yielding 

exp[ 79 (2iV+l) d ] (iV^+oo), (18) 

which completes the proof of Lemma 2. □ 

It follows from (12) and (18) that there is a transition from a regime where 
\im.N^ +OQ (£N(x, g)) < +oo to a regime where \im.N^ +00 (£N(x, g)) = +oo as g increases 
past g — 1. The value g — 1 (or, more exactly, A = 1/t) is the counterpart of what we called 
the "critical coupling" in [5]. 

Note also that (£]y(x,g)) does not depend on x. In the following we will write (^(fl 1 )) 
for (£ N (x,g)). 



(£ N (x,g)) 



2vr 



(2N + if 2 



1+0 



J- 

\ N d / 2 



IV. INTERMITTENCY OF € N (x,g) AND LOSS OF ERGODICITY 



In this section we investigate the intermittency of S^(x,g) as N — > +oo, assuming er- 
godicity for g < 1 where liniAr^ +00 (£Ar(g)) exists. 



Proposition 1 IfV g <1, one has 

Jim / 5jv(a;, 5) d d x = Jim (£^(3)) = -J_, 
almost surely, then M g > 1, 



and 



almost surely. 



N 



lim / £ N (x, g) d x = +00, 
n^+oo J A 



lim 1 / S N (x,g) d d x = 0, 
■-+00 (£jvU0) A 



(19) 



(20) 



(21) 



[For |A| 7^ 1, the left-hand side of Eqs. (19) to (21) is divided by |A|]. Assuming ergodicity 
for g < 1 [Eqs. (19)], one finds that f A £]y(x,g) d d x is not asymptotic to (S^(g)) as A" — > +00 
if g > 1 [Eq.(21)]. In this asymptotic sense one can say that ergodicity is lost as g increases 
past the value g — 1. Note also that (19) ensures that the assumption we made to prove (7) 
is fulfilled: use E^^x^g) 6 = £^(x,ge) and take e < 1/g. 



Proof. From (10) it follows that for every N > and every realization of the 9 n , 
J* a £at(x, (?) d d x is a non decreasing function of g. Thus, for every g > 1 and e > 0, 

/ £ N (x, g) d d x > / £ N (x, 1 - e) 

A </A 

and by Eq. (19), 

liminf /" ^(x, #) d d x > lim / £jv(x, 1 - e) d d x = -, (22) 

Ar - + +°° A A->+oo J A £ 

almost surely. Letting e yields (20). 

We now prove the limit (21). From the control (A. 6) with a = 3/4 it follows that for N 
large enough, 



J A 

< exp [2#(2iV+ l) d/2 ] , (23) 



5fsup |SW(aO 



almost surely, and by Lemma 2, 

v 111 ? / 9) d " x = °> ( 24 ) 

a->+oo {£ N (g)) J A 

with probability one, which completes the proof of Proposition 1. □ 

Proposition 2 (transition to intermittency) Under the same ergodicity assumption as 
in Proposition 1, g th = 1. 

Proof. First, we prove (? tft > 1. If £ N (x,g) is intermittent, then for every integer p > 1 
and almost all the realizations of £^(x,g), there exists / p _i(jV), with / p _ 1 (iV) — > +oo as 
iV — > +oo, such that 

/ £a^,#) p ^x~ / ^ Af (x,^) p l|^ Jv(;ri9) | >/p _ l(iV) rf d x, (25) 

</ A JA 

as iV — > +oo. From Holder's inequality and 

^ le N {x,g)>MN)d?x < — !^ £ N (x, g) d d x, 
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one gets, Ve > 0, 



£ N (x, g)le N (x,g)>f (N) d d x 



< 



£ N (x,g) 1+e d d x 



l+s 



d d x 



< 



£ N (x,g(l + e)) d d x 
1 



A 

i 

l+e 



l+e 



fo(N) 



l+e 



l£ N (x,g)>fo(N) 

/ ls N (x,g)>MN) d d X 
J A 

£ N (x,g(l + e)) d d x 



l+e 



(26) 



l+e 




l+e 




/ £ N (x,g)d d x 






J A 





Take g < 1/(1 + e). By (19), the two brackets on the right-hand side of (26) are almost 
surely bounded and J A £n(x, g)le N (x,g)>f (N) d d x — > almost surely as iV — > +oo. This is in 
contradiction with (25) for p = 1. Thus, for every g < 1/(1 + e), £N(x,g) is not intermittent 
and by taking e > arbitrarily small one obtains g t h > 1- 

We now prove ^ < 1. To this end we prove that (2) is fulfilled if g > 1. Using Holder's 
and Jensen's inequalities successively, one finds that for every integer p > 1 and V0 < e < 1, 



/ £ N (x,g) p d d x 

J A 

S N {x,gf v+1 ^-^d d x 
[ S N {x,g) {p+1){l - £) d d x 

J A 



1 




p+l 








1 




p+1 


1 







p(p+l) 



p(p+l) 



(27) 



< 



£ N (x,g) p+1 d d x 



p+l 


'J, 







p(p+l) 



s N (x, g y^s N (x, g y £ d d x 
s N {x, g y( p+ V £ d d x 

£N(x,gp(p+ l)e) d d x 
£ N (x,gp(p + l)e) d d x 

Ua J Ua 

which gives, 

i _i 
r \ p+i \ p p 

£ N (x,g) p d x 

£N(x,gp(p+l)e) d d x 
£ N (x,gp(p+l)e) d d x 

For every g > 1, take < e < [gp{p + l)] -1 - By (19) and (20), one has, with probability 

(29) 



> 



p+i 



£ N (x,g) p+1 d d x 
£ N (x,g) p+1 d d x 
£ N {x,g{p + l))d d x 



A 

£ 

P+l 



"p(p+l) 



(28) 



"p(p+i) 



one 



N^+oo J A 
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and 

lim / £n(x, q(p + 1)) d d x = +00. 
Injecting (29) and (30) into the right-hand side of (28), one gets 



(30) 



£ N (x,g) p d d x 



£ N (x, g) p+1 d d x 



(31) 



almost surely. It remains to prove the first inequality (2), which is immediate by (20). 
Thus, (2) is fulfilled and £^{x,g) is intermittent for every g > 1. This implies g t h < 1, 
which completes the proof of Proposition 2. □ 



V. SUMMARY AND PERSPECTIVES 



We have investigated intermittency of SN(x,g) = exp[g\SN(x)\ 2 ] as N — > +00, when 
Sn(x) is given by the Steinhaus series (9). Assuming ergodicity of £j^(x,g) as N — > +00 
for g < 1, where limjv^+oo^AKfl 1 )) exists, we have proved the existence of a transition to 
intermittency as g increases past the threshold g t h = 1 (Proposition 2). This transition goes 
together with a loss of ergodicity in the sense of a transition from (assumed) ergodicity at g < 
g th to a regime where Xm\ N ^+J^k\{£ N (g))Y x § K £ N {x,g) d d x = at g > g th (Propositions 
1). Proving ergodicity of £js[{x,g) as iV — > +00 for </ < 1 is another problem that we are 
unable to solve at the present time. 

The next step toward a study of the solution to (1) will consist in allowing for a time 
dependent Sn in the simpler Laplacian free case (\m\ = +00). What is to be expected in this 
setting can be conjectured in view of the results obtained in this paper. Indeed, it is easily 
seen that the intermittency threshold of Proposition 2 corresponds to the critical coupling 
defined as the smallest g at which 0)) would diverge if Sn(x) was a Gaussian r.v. 

with (S N (x)) = (S N (x) 2 ) = and (\S N (x)\ 2 ) = 1. Intuitively, this could have been expected 
from the CLT according to which, for any fixed x G M. d , (9) tends in law to such a Gaussian 
r.v. as N — > +00. Now, if Vx G R d , Sn(x, t) tends in law to a Gaussian random function 
of t as iV — > +00, it is not unreasonable to conjecture that the intermittency threshold for 
£isr(x,t) should be given by the critical coupling in this case too. Namely, for any given 
t > 0, 

Xth{t) = iwr 
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where /J,i(t) is the largest eigenvalue of the covariance of Sn(x,t) for < r < t. Proving 
this conjecture goes through the resolution of specific technical problems inherent in £n{x, t) 
being now a functional of S]y(x,t). This will be the subject of a forthcoming paper. 

We will then be ready to tackle the study of intermittency of the solution to (1) with a 
finite m. This more difficult problem will presumably require the use of the distributional 
formulation of [5]. 

APPENDIX: CONTROLLING THE EXCURSION OF \S N {x)\ 

Pave A with rf-dimensional cubes, A,, of length t < 1, with t~ x G N and 1 < % < £~ d . Let 
ij G A denote the center of Aj. For every x G A, one gets 



\S N (x)\ 2 -\S N ( Xi 



< 



< 



E 

2 



(2N+1) 



— e 



^2in(n—m)-x 2in(n— m)-Xi 



(2N+1) 
(2N+l) d 



— ^ \sm[7r(n - m) ■ (x - Xi 



E 



\n — m\\x — Xi 



Since \n — m\ < 2N\fd and \x — Xi\ < £\/d one has Vi e Aj, 

\\S N (x)\ 2 - \S N ( Xi )\ 2 \ <4ndW(2N + l) d . 
Fix 1/2 < a < 1 and take 

r 1 = int(47rd + 1)N{2N + l) 2 ^ 1 -"). 

Then, (A.l) is bounded by 

\\S N (x)\ 2 -\S N (x i )\ 2 \<(2N + l) d ^- 1 \ 

Let zgC and write t = \z\ and if = Arg(z) — Arg(5 , 7v(^))- For every x G A, one has 

2* 



(A.l) 



(A.2) 



(A.3) 



(exp[z<Sjv(a;)* + c ■ c-]) = exp 



(2N+l) d \n I c 



(2JV + l) d A 



and 



(exp[zSV(^)* + c • c-]) = (exp[2t\S N (x) | cos </?]) 

> (exp[2t|Sjv(aOI cosy9]l_ V 3<^ /3 l|^ W |> (2W+1) d( Q -i/2)) 

> iexp [t(2iV + l) d (*-V2)] p (1^(^)1 > (2AT + l)^ 1 / 2 )) . 
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Thus, 



P(\S N (x)\ > (2N + lY^-W) 



< 3 exp 



(2iV + l) d ln I c 



2t 



(2N + \y/ 2 



t(2N + l) d («-V2) 



<3exp [t 2 - t(2N + l)^ 1 ^] , 



where we have used the inequality In Io(2s) < s 2 , and by taking t = (2N + l) d ( a 1 / 2 ) /2 one 
gets 



--(2AT+ i)^- 1 ) 



(A.4) 



^ (\Sn(x)\ > (2N + l) d ( Q -V2)) < 3exp 

From (A.3), it follows that if 1^(^)1 < (2iV + l)^" 1 ^) for every % < t d , then |SW(aO| < 
y/2(2N + lyi*- 1 / 2 ) for every x E A. Therefore, using (A.4), 



> (27V +1) 

i — 1 



d(a-\/2) 



} 



< 3r d exp 



■-(2N +l) d{2a - L) 



(A.5) 



Now, by (A. 2) £ d diverges algebraically in A" as A" — > +oo, and since a > 1/2 the right- 
hand side of the last inequality (A.5) tends to zero faster than any power of A" as A" — > +oo. 
Thus, 

V P (sup ^(aOI > ^(2^ + l)«*(«-i/2) j < +00) 



and by the Borel-Cantelli lemma, 



P Mini sup sup |SW(a;)| < ^2{2N + l)^" 1 / 

\ A^+oo a;eA 



/ 2 ) = 1. 



(A.6) 
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